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<N . Abstract 

r~| . We consider the case of a pair of particles initially in a superposition state corresponding to a 

separated pair of wave packets. We calculate exactly the time development of this non-Gaussian 
C^ ' state due to interaction with an arbitrary heat bath. We find that coherence decays continuously, 

CT^. as expected. We then investigate entanglement and find that at a finite time the system becomes 

separable (not entangled). Thus, we see that entanglement sudden death is also prevalent in 

> 

f^ ' continuous variable systems which should raise concern for the designers of entangled systems. 
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For continuous variable systems "entanglement sudden death" [1], that is, complete ter- 
mination of entanglement after a finite time interval, has been demonstrated for the special 
case of a pair of particles in a Gaussian state [2]. Those authors use a master equation and 
the necessary and sufficient criterion for separability of such states developed by Duan et 
al. [3]. Here, we present a more general model by considering the case of a widely separated 
pair of particles initially in a superposition state corresponding to a displaced pair of wave 
packets. We use a method that allows us to calculate exactly the time development due to 
interaction with an arbitrary linear passive heat bath [4] . We find first of all that coherence, 
defined as the relative amplitude of the interference pattern, decays continuously but very 
rapidly. Next we consider entanglement and find that after a finite time the system becomes 
separable, showing that "sudden death" of entanglement occurs for this system as well. 

The method is based on the general prescription described in an earlier publication [4] in 
which a system is put in an initial state by a measurement applied to the equilibrium state 
and after a finite time is sampled by a second measurement. A key formula is the expression 
for the Wigner characteristic function given in Eq. (6.5) of [4]. For a two particle system, 
this formula takes the form: 

W{Q„ Pi; g^, P2) = ^^^^ {Piljm ' ^^^ 

where the initial measurement is described by 

/(l) = /(a;i(0)-xi,X2(0)-X2), (2) 

in which f{xi,X2) is the c-number function describing the initial measurement while Xi{t) 
and X2{t) are the time-dependent Heisenberg operators corresponding to the displacement 
of either particle: 

a:,•(t) = e'^*/'^a;,(0)e-'^*/^ (3) 

Finally, in this formula the brackets indicate expectation with respect to the state of the 
system in equilibrium at temperature T, 

Tr < OekT I 

(O) = \ ,/ . (4) 

Tr I efcT I 

Here we emphasize that in Eqs. (3) and (4) H is the Hamiltonian operator for the system 
of the pair of particles interacting with the heat bath. 



In order to evaluate this formula we make the key assumptions that particles are linear 
oscillators coupled to a linear passive heat bath and that within the bath the particles are 
widely separated so that we may ignore bath-induced interactions. These assumptions imply 
that Xi{t) and X2{t) independently undergo quantum Brownian motion. We can now repeat 
the discussion leading to Eq. (6.43) of our earlier publication, [4] to obtain 

/ J!U]\p~i{^i{t)Pl+mxi{t)Qi+X2{t)P2+mx2{t)Q2)/h f(l)\ 
n=l 

X / dx[ J dx'j\x[ + y, 4 + y)/K - y> 4 - y) 



where (x^) and (i^) are the mean squares of the displacement and velocity, the same for 
either particle, and we have introduced 

^" = — ^7"^\ ^' ^" = '^^" + '^^Qn- (6) 

[X ) 

Here G = G{t) is the Green function and c = c(t) = ^{x{t)x{0) +x{0)x(t)) is the correlation 
function, again the same for either particle. 

These expressions are valid for any measurement function. We now specialize to the case 
of a pair of particles initially in a superposition state corresponding to a separated pair of 
wave packets, with measurement function of the form: 



f{Xl,X2] 



1 



V47ra2(l + e-'^'/4-^) 



exp 



{xi - d/2f + {X2 + d/2) 



4a2 



+ exp 



(xi + d/2Y + (X2 - d/2) 



(7) 



4^2 

Here we emphasize that the wave packet separation d is arbitrary and should not be confused 
with the separation of the particles in the bath, which is large. 

With this measurement function the integrals in the expression (5) are standard Gaussian 



[4]. Putting the result in the expression (1) for the Wigner characteristic function we find 

W{Qi,Pi;Q2,P2;t) 



^ {^') {P^ - 1^) + m' {i') Ql + ^Ll 



n=l 

X S K ^^ -■ (8) 

i + eXp <j^ 4a2((x2>+a2) J 

where in order to center the state at the origin we have put aji = 0:2 = 0. 

This expression becomes simpler in the free particle limit : {x^) — ?> 00. In this limit 

W{Qi,Pi-Q2,P2;t) 



exp 



An(Pf + P|) + 2A^2{QiPi + Q2P2) + A22{Ql + Ql) 



2h 



2 



cos (^^#1^ + exp {-i^l cosh [^(^^-^2)+^c(Q.-Q.)]. 

X ^ r r^^ -^ (9) 

in which we have introduced 

2 n^G' 



All = (r^ + s + 
A12 



4a2 ' 
ms fi^mGG 



2 4(t2 ' 

^. = m=(i^> + ^^. (10) 

In these expressions s = 2 {{x^) — c) = {{x(t) — x{0)) ) is the mean square displacement 
and as above G is the Green function. 

The Wigner function is the inverse Fourier transform of the Wigner characteristic func- 
tion: 

1 



W{qi,pi;q2,P2;t) 



2(1 - e-'^V4-2 



Woiqi - -,Pi; t)Wo{q2 + -,P2; t) 



+Wo{qi + -,Pi;t)Wo{q2 - ^,P2;t) 
+2e-^^'^Woiqi,Pi;t)Woiq2,P2;t)cos^qi-q2,Pi-p2;t)]. (11 



Here Wq is the Wigner function for a single particle wave packet, 

uif ,\ 1 / A22q^ - 2Ai2qp + Aup^ 

"^'^'^'-''^ ^ 2nVAiiA.2-Air'V 2iAiiA22 -Al,) ^ ^''^ 
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while the phase $ is given by 



"-^ — 




A11A22 - A\^ 


4a2 


(^11- 


- 4^2 K^22 


4^2 ) 1^12 


rP-mGG\2 j2 
- 4a2 J " 



^, ^. {GA22 - mGAi2)q + (mGAn - GAi2)p M 
and the quantity A by 

^ ^ " AnA22 - Al, 4^- ^^^^ 

We note that each of the first two terms in brackets in the expression (8) for the Wigner 
function corresponds to the product of independently propagating packets. We call these the 
direct terms. The third term is an interference term. We emphasize that we have assumed 
that the particles are widely separated within the bath so there is no coupling between 
them. The presence of this interference term is therefore a purely quantum mechanical 
phenomenon. 

The Wigner function is a quasiprobability distribution, not directly observable. A phys- 
ical observable is the probability distribution, obtained by integrating over the momentum 
variables: 



1 

P[qi,q2]t) 



2(1 - e-'^V4a 



2^ 



Po(gi--,t)Po(g2 + -,t) 



+Po(gi + 2>^)^o(g2--,t) 

+2a(t) exp {-^} Po(gi; t)Po(g2; t) cos I """^If^-/'^ 



(15) 



Again, the first two terms are direct terms corresponding to independently propagating wave 
packets with 

Po(g;t) = ^^=exp( ^ 1 , (16) 

the probability distribution for a single wave packet centered at the origin. The third term 
is an interference term. Viewed in the (g'1,5'2) plane, the direct terms are seen as a pair 
of peaks centered at {qi,q2) = (f, — f) and {qi,q2) = ("f' 2) ^^^ spreading in time as the 
width All increases. The interference term is seen as a spreading peak centered at the origin 
and modulated by the cosine term. The quantity a(t) is the ratio of the geometric mean of 
the direct term to the factor multiplying the cosine in the interference term and is therefore 
a measure of the visibility of the interference. We find 



This quantity is initially unity and, for d large, diminishes rapidly to a very small asymptotic 
value. This is the familiar phenomenon of decoherence of a superposition state. But never- 
theless interference is present for all times, albeit with a small amplitude. Our point here is 
that there is no sudden death of coherence as indicated by the presence of the interference 
term. 

We turn now to the question of entanglement. A two-particle state described by a density 
matrix p is said to be separable (not entangled) if and only if p can expressed in the form 

P = J2P^PA^)Pji2), (18) 

j 

in which Pj{l) and Pj{2) are projection operators into states of particles 1 and 2, respectively, 
and the pj are positive. In our case we seek to express the density matrix elements in the 
form 

{x[,X2\p\xi,X2) = / (fai / (fa2P {ai,a2) (J)ai{x'i)(l)l^{xi)(f)a2{x2)(j)*a2{x2), (19) 
where the 0's are what we might call strong form coherent wave functions: 

(j)a{x) = [2TTaQ) exp<^ 4^ (^ ^ ^) + ^ ~ ^2^ f ' ^^^^ 

with the state labelled with the complex number 

l-i6o_ , .ao_ ,2 dxdp 

a = x + t--p, da = — -— . (21) 

2cto ^ 2h ^ ^ 

This is clearly of the form (18) with the sum replaced by an integral, so if this expansion 
exists and P{ai,a2) is everywhere positve the state is separable. The expression (19) is 
reminiscent of the Glauber-Sudarshan P-representation, [6] but in that representation the 
0's are coherent states, which are expressed in terms of the ground state of an oscillator, 
or equivalently a minimum uncertainty state, [5] shifted in position and momentum. If in 
the wave function (20) we set the parameter 6q equal to zero we have such a coherent wave 
function. On the other hand, if 60 is not zero, the wave function (20) minimizes the strong 
form of the uncertainty relation: [7, 8] 

{{x-xf) {{p-pf) - /i--^)iP-P) + iP'P)i--^)V > ^. (22) 



2 / ~ 4 

It is not difficult to show that wave function (20) satisfies this as an equality. 
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Next, we recall the relation between the Wigner characteristic function and the density 
function matrix elements: 



W{Qu Pi; Q2, P2) = J dqij rfgse-^^'^^^^+^^^^^Z'^ ^^i - ^, ^2 - ^ |p| gi + ^, ^2 + ^ ) . 

(23) 
Using the expansion (19) of the density matrix elements, this becomes 

W{Qi, Pi;Q2, P2) = J d^ai J Sa^P (ai, a^) J dqi<j)^, (gi - ^)0;^ (gi + ^)e-*'?^^i/^ 

X jdq,<P^Aq, - ^)0:,(g2 + ^)e-^'^^^^/^ (24) 

With the explicit form (20) of the coherent state, we see that 

,, _ |)«, 4)e--/« . e-<-*- exp {-f! - M?^ - (i±f£} ^ (25) 

Therefore, the expression (24) can be written 



i=i,2 



-^ —, '-^ '- } W{Q^, Pv, Q2, P2). (26) 



This is just the Fourier transform of the P-function, which will exist if the inverse transform 
exists. From an inspection of the Wigner characteristic function (9) for our superposition 
state, we see that convergence of the inverse transform will be dominated by the exponential 
factors and will therefore exist if the quadratic form 

is positive definite. Since the parameters Gq and b^ are arbitrary we can first choose 8q 
to make this quadratic form diagonal and then choose o"o to minimize the product of the 
diagonal elements. The corresponding optimum values are 



('^o)opt-^-, K)op,-y ^^ • (28) 

With this choice we find for the diagonal elements that the diagonal elements of the quadratic 



form (27) are given by 



^1 = \ -r- \/^ii^22 - W^?2 + -r 



A 



22 



^2 = \/^(^^4;^-y^?2 + ^ ■ (29) 

It is not difficult to see that these are positive at all times. Thus the expansion (19) exists 
at all times. 

Next we consider the positivity of P{ai,a2) ■ With the opimum values (28) of the 
parameters in (26) we form the inverse Fourier transform. The integrals are all standard 
Gaussian and the result can be written in the form 



P{ai,a2) 



n exp <^ ^^^^ ^^^^ ^^1 ^^,^^^ ^^^^ J -^ 

7r2ini22(l + exp{-^}) 

, , ^ n^G^ h^a'm^G^ a^ \ £ ] , {x, - X2) d 

exp < 1 = ^ — — ■ > cosh 

^' 4(TMn 4^22 AiiJ^a^j 2A, 

(hGd {X1-X2) , hmGd{pi-p2) , ,„^. 

+ cos = \ = . (30) 

\^ 4^2^11 4^2^22 ' 

The ffist line in Eq. (30) is a positive factor, so P {ai, 0:2) is positive if the remaining factor 
is positive. Clearly this will be the case for all {xi,pi,X2,P2) if and only if 

G{t) = 1 . . ^ > 0. (31 

4^2^11 4^2^22 An 

At short times G{t) = t/m and s{t) = {x^) t^. With the expressions (10) for Au^ A12 and 
^22 and these in turn in the expressions (29) for An and A22 we find 

CO, ^ -^ . _i±li±i^!Z^, (32) 



where A = fijm^lTW) is the deBroglie wavelength. Not surprisingly C(0) is always negative, 
since the initial state is formed with a projection operator (7) corresponding to a necessarily 
entangled state. 

At very long times, the behavior of G{t) and s(t) depends upon the bath parameters. [9] 



As an illustration we consider the Olimic model for which at long times 

G{t) ~ C\ 

s{t)r^ ^\og^, T = 0, 
irC, rn 

2kT 
s{t) —t, T > 0, (33) 

where ( is the Ohmic friction constant. With this it is easy to see that for this Ohmic case 
at long times C{t) ~ 1. Clearly there must be an intermediate time at which C{t) changes 
sign and the state becomes separable. For example, in Fig. 1 we plot C{t) versus '-ft for 
the single relaxation time (r) model [4] at zero temperature, where r = '~f^^/Q and 7 is the 
Ohmic relaxation time. There we see that the change of sign occurs at '~ft ^ 6. In general, 
most other bath models (that is, models with colored noise [9]) show similar behavior. 
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FIG. 1: C{t) versus 'jt for the single relaxation time (r) model at zero temperature, where r 
7^^/6 and 7 is the Ohmic relaxation time. We note that separability occurs for 'yt ^ 6. 
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